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Abstract 

We investigate the matching of continuous gravitational wave (CGW) signals in 
an all sky search with reference to Earth based laser interferometric detectors. We 
consider the source location as the parameters of the signal manifold and templates 
corresponding to different source locations. It has been found that the matching of 
signals from locations in the sky that differ in their co-latitude and longitude by 7r 
radians decreases with source frequency. We have also made an analysis with the 
other parameters affecting the symmetries. We observe that it may not be relevant to 
take care of the symmetries in the sky locations for the search of CGW from the output 
of LIGO-I, GEO600 and TAMA detectors. 
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1 Introduction 

The first generation of kilometer-scale gravitational wave (GW) laser interferometric de- 
tectors with sensitivity in the frequency band 10 Hz to few kHz and ultra cryogenic bar 
detectors sensitive at frequencies around 1 kHz will start collecting data soon. The TAMA 
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300 (Tsubona 1995) has already done the first large scale data acquisition (Tagoshi et al. 
2001), while LIGO (Abramovici et al. 1992) and GEO600 (Danzmann 1995) have recently 
carried out their first science observations. VIRGO (Bradaschia et al. 1991) may become 
operational in couple of years. Also, an eighty meter research interferometer ACIGA (Mc- 
Cleland et al. 2000) near Perth, Australia is under construction, hoping that it may be 
possible to extend it to multi-kilometer scale in the future. 

At present, majority of searches are focussed in the detection of chirp and burst signals. 
However interest for the search of continuous gravitational wave (CGW) signals from the 
output of detectors is growing (Jaranwoski, et. al 1998, Jaranowski & Krolak 1999, 2000; 
Astone et al. 2002; Brady et. al. 1998, Brady & Creighton 2000) due to the possibility 
of enhancing the signal-to-noise ratio (S/N) by square root of observation time y/T obs . An 
optimistic estimate suggests that the Earth based laser interferometric detectors may detect 
such signals with an observation time of 1-yr. 

The strength of the CGW largely depend on the degree of long-lived asymmetry in the 
source. There are several mechanism for producing such an asymmetry (Pandharipande et 
al. 1976; Bonazzola & Gourgoulhon 1996; Zimmermann & Szedenits 1979; Zimmermann 
1980). The estimates of the asymmetry in neutron stars shows that the amplitude of CGW 
may be < 10 -25 . Hence, long integration will be required to get the signature of signal. But 
this in turn induces several other problems viz. Doppler modulation and non-stationarity 
of the detector noise. Consequently data analysis becomes more harder. However, Doppler 
modulation will provide the information of position of the source in the sky. 

It has been realized that with present computing power, coherent all- sky full frequency 
search in the bandwidth (BW) of the detector of a month long data is computationally 
prohibitive. Some computationally efficient alternative approaches have been suggested 
viz. tracking and stacking (Brady and Creighton 2000; Schutz 1998). Tracking involves the 
tracking of lines in the time-frequency plane built from the Fourier transform (FT) of one 
day long stretches of data while stacking involves dividing the data into day long stretches, 
searching each stretch for signals, and enhancing the detectability by incoherently summing 
the FT of data stretches. Also, accurate modeling of GW form, optimal data processing and 
efficient programming are an integral part of all sky-search. 

The basic method to analyze the detector output to get the signature of GW signals rely 



2 



on how efficiently one can Fourier analyze the data. Fourier analysis of the data has the 
advantage of incorporating the interferometers noise spectral density. The problem for the 
search of CGW largely depend on how accurately one can take into account the translatory 
motion of the detector acquired from the motions of the Earth in Solar system barycentre 
(SSB) frame. It has been shown (Srivastava and Sahay 2002a,b) that amplitude modulation 
will only redistribute the power of the frequency modulated (FM) signal in five frequency 
bands / ± 2f rot , f, f ± f rot , where / and f rot are the frequencies of the FM signal and the 
rotational frequency of the Earth respectively. Hence it is sufficient to consider only FM 
signal for the analysis of the matching of signals from different locations in the sky . 

The study of the matching of signals from locations in the sky that differ in their co- 
latitude and longitude by it radians has been made for 1-d data set (Srivastava and Sahay 
2002c). They observed symmetries in the sky locations. In this paper we report that the 
observed symmetries are not generic. We study this problem for longer observation time 
and the other parameters which affect the symmetries. In the next Section we briefly review 
the FT of FM signal. In Section 3, using the concept of fitting factor (FF), we investigate 
the matching of signals in an all sky-search with reference to the Earth based laser inter- 
ferometric detector by considering the source location as the parameter of signal manifold 
and templates corresponding to different source locations for longer observation time. We 
present our conclusions in Section 4. 

2 Fourier transform of the frequency modulated continu- 
ous gravitational wave 

The time dependence of the phase of a monochromatic CGW signal of frequency f Q ob- 
served at detector location is given as (Srivastava and Sahay 2002a) 
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where 

V = 2nf — sin a(cos/3 (sin# cos e sin + cos 6* sin e) — sin/^ sin 9 cos0) , 

Q = 2nf ^ sin a(sin / 9 (sin6 ) cose sin0 + cos 6* sin e) + cos/9 sin 9 cos0) , 

A/" = VP 2 + Q 2 , 
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where R e , R se , w rot and w orb represent respectively the Earth's radius, the average distance 
between the centre of Earth from the origin of SSB frame, the rotational and the orbital 
angular velocity of the Earth, e and c represent the obliquity of the ecliptic and the velocity 
of light, a is the colatitude of the detector. Here t represents the time in s elapsed from 
the instant the Sun is at the Vernal Equinox and @ is the local sidereal time at that instant, 
expressed in radians. 9 and <fi denote the celestial colatitude and celestial longitude of the 
source. These coordinates are related to the right ascension, a and the declination, 5 of the 
source via 

cos 9 = sin 5 cos e — cos 5 sin e sin a 



sin 9 cos <fi = cos 5 cos a 

sin 9 sin = sin 5 sin e + cos 8 cos e sin a 
The two polarisation states of the signal can be taken as 

h+(t) = h 0+ cos[$(t)] 



(4) 
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h x (t)=h 0x sin[$(t)] 
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(6) 



h Q+ , h Qx are the time independent amplitude of h + (t), and h x (£) respectively. 



For the analysis of the matching of signals from different locations in the sky, it is 
sufficient to consider either of the polarisation states of FM signal. We consider the + 
polarisation of the signal of unit amplitude. Therefore, the FT for a data of T obs observation 
time may be expressed as (Srivastava and Sahay 2002b) 
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J stands for the Bessel function of first kind. The computational strain to compute h(f) 
from equation © can be reduced by »s 50% by using the symmetrical property of the 
Bessel functions, given as 
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X = sm(JZ + Q-m-n/2) 

y = cos(ft + Q - mvr/2) ^ 

U = sin cos m(£ — 5) — ™ {cos z/£„ sin — 5) — sin m5} 

V = cos cos m(£ D — 5) + — sin i/£ D sin m(£ Q — 5) — cos m<5 

Equation © contains double infinite series of Bessel function. However, we know that 
the value of Bessel function decreases rapidly as its order exceeds the argument. From 
equations ©, it can be shown that the order of Bessel function required to compute h(f) 
in the infinite series are 



3 Matching of the continuous gravitational wave 

Matched filtering is the most suitable technique for the detection of signals whose wave 
form is known. The wave forms are used to construct a bank of templates, which represent 
the expected signal wave form with all possible ranges of its parameters (source location, 
ellipticity, etc.). In an all sky search by match filtering, Srivastava and Sahay 2002c ob- 
served symmetries under the following transformations 

9 T ► 7T - 6 T 0<9 t <ti, (13) 

(j) T > 7T — (f) T < < 7T , (14) 

<pT ► 37T — (f)T 7T < <pT < 27T. (15) 
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The above observations are made by studying few cases for 1-d data set. Hence, it will 
be important to understand and check the generic nature of the matching of signals under 
the above transformations. To study this problem we use the formula for FF (Apostolatos 
1995) which quantitatively describes the closeness of two signals, given as 

PP = (h(f)\h T (f;e T ,<p T )) 

y/(h T (f;9 T ^ T )\h T (f;9 T ,<p T ))(h(f)\h(f)) 

where h(f) and hr(f; 9 T , <p T ) represent respectively the FTs of the actual signal wave form 
and the templates. The inner product of two waveform hi and h 2 is defined as 

= 2 r mMfp^MnMfi df 

Jo &n\j) 



oo 



= 4 / m (17) 

Jo b n {J ) 

where * denotes complex conjugation, ^denotes the FT of the quantity underneath (a(f) = 
J°° a(t)exp(— 2nift)dt) and S n (f) is the spectral density of the detector noise. In our 
analysis we assume the noise to be stationary and Gaussian. It is remarked that to compute 
the inner product of the signals one would require the BW of the Doppler modulated signal. 
From the Doppler shift (Srivastava and Sahay 2002a), the BW of Doppler modulated signal 
may be given as 

BW « (1.99115 x 10~ 4 sin 9 + 3.09672 x 10' 6 )f o . (18) 



3.1 Celestial colatitude 

Let us consider that GEO600 detector (the position and orientation of the detectors can 
be found in Jaranowski et al. 1998) receive a CGW signal of frequency f a = 0.1 Hz 
(unreasonably low frequency has been chosen for illustrative purposes limited by available 
computational power) from a source located at (8, 4>) = (25°, 20°). In order to evaluate the 
matching of the signals in colatitude, we first maximize FF over by choosing = <\> T = 
20°. Now, we wish to check the symmetries in colatitude represented by equation (fT3t for 
the data set T obs = 120 d. For the purpose we maximize the FF over 9 by varying Q T in 
discrete steps over entire range i.e. 0° to 180°. For the present case it is sufficient to take 
the ranges of k and m as 1 to 345 and —3 to 3 respectively and BW = 20.1954 x 10~ 6 Hz. 
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The results so obtained are shown in Figure ©. To establish the observed symmetries, 
we similarly compute the FF for T obs = 1, 2, 3. ...365 d and observe that the matching of 
signals remains almost same. However, due to obliquity of the ecliptic, the variation in the 
matching of signals will depend on source frequency, colatitude and detector position and 
orientation. We check the dependence of the FF on this parameters. The result so obtained 
for different Earth based laser interferometric detectors are shown in the Tables ©, & © 
and Figures ©& @. 

The analysis of the matching of signals in 9 space shows that 

(i) for fixed f a , the FF 

(a) is independent of T obs and 4>. 

(b) does not vary significantly with the variation of source location, detector posi- 
tion and orientation [Tables ©, and ©]. 

(ii) the FF falls with the source frequency [Figures © and ©]. From the figure we find 
that it may not be relevant to take care of the symmetries in the sky locations for the 
search of CGW from the output of LIGO-I, GEO600 and TAMA detectors whose 
lower cut off frequency is 40/75 Hz (Owen and Sathyaprakash 1999). 

(iii) the approximate fall of FF based on the figures ©and © may be given as 



where A a , A\, A%, A 3 are constants given in Table ©. 

3.2 Celestial longitude 

The Doppler shift due to the motions of the Earth is mainly depend on the colatitude and 
source frequency and have very less dependence on the longitude. Consequently, grid 
spacing of the templates for matched filtering for an all sky search will insignificantly de- 
pend on longitude (Brady & Creighton 2000). Keeping this in view we similarly check the 
matching of signals under the transformation given by the equation (fT4l) . We chosen the 
LIGO detector located at Livingston, selected a data set T obs = 120 d, (9, cp) = (0.5°, 40°), 
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Figure 2: Fall of FF with f . 
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Figure 3: Fall of FF with f a . 
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Table 1 : Matching of the signals of frequency 1 Hz under the transformation represented 
by equation (fT3l for GEO600 detector. 
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Detector 


A 


Ai 


A 2 


^3 




xlO" 3 


xl0~ 5 


xlO" 5 


xiir 7 


GEO600 f/3 - °) 


1000.02 


124.524 


129.030 


137.774 


GEO600 (/3 _9o^ 


997.796 


175.426 


248.072 


282.478 


LIGO Hanford^ _ <^ 


998.450 


225.359 


172.731 


221.187 


LIGO Hanfordfa-gno-) 


996.668 


266.811 


324.073 


440.729 


VIRGO (/3o=0 o) 


998.548 


243.480 


191.503 


259.026 


VIRGO {/3o=900 ) 


996.064 


316.288 


362.630 


529.531 


TAMMA300 (/3o= oo) 


997.914 


342.746 


249.527 


390.874 


TAMMA300( / g o =90«») 


995.744 


371.193 


458.856 


756.536 


LIGO Livingston ( - /3o=0o - ) 


998.815 


396.783 


285.368 


478.932 


LIGO Livingstori( /3o=90 o) 


995.622 


399.087 


516.416 


904.211 



Table 2: Coefficients of the fall of FF with f a under the transformation represented by 
equation O for (3 = 0° and 90°. 

f = 5 Hz. In order to compute the FF, we first maximize equation (fToT) over 9 by se- 
lecting 9 = 9t = 0.5°, followed by maximization over <p in discrete steps over its entire 
range, 0° < <fr < 360°. The result so obtained is shown in Figure ©. We also check the 
mismatch of the signals for different 9, <p and /„ by computing the FF for the data set of 
T b s = 1,2, ....25/100 d. The results so obtained are shown in Figures ©, ©, and © 
respectively. Almost same behavior has been observe for the transformation represented 
by the equation (fT31) . 

From the figures we note that the matching of the signals in longitude decreases with 
T bs, fo, 9 and </>. However, behavior of the matching of signals are similar in nature and 
may be represented by 

FF = B - B 1 T obs + B 2 T 2 obs - B 3 T* bs + B 4 T^ bs ; T obs = 1, 2, ....25/100 d (20) 

Where, B a , B x , B 2 , B s , B 4 are the constants given in Table ©. The equations ( fT9l 
and (|2H does not represents the oscillatory part of the figures. This oscillatory behavior is 
more or less typical in waveform that match well or bad depending on their parameters and 
is low enough in comparison to the threshold of the detection of GW that one may choose. 
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19567.0 
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1032.67 


13220.3 
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56298.5 


54579.0 
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15 






1034.40 


16673.7 


21157.0 


103289.0 


122385.0 


1 






1/181 


179/359 


1031.20 


2896.18 


722.052 


671.331 


149.659 








50/230 


130/310 


1028.29 


2216.41 


452.688 


341.054 


61.1078 




0.5 


0.5 


70/250 


110/290 


1024.32 


1473.47 


225.270 


125.891 


16.4451 








80/260 


100/280 


1019.87 


906.855 


101.752 


41.5159 


3.84672 








85/265 


95/275 


1012.58 


449.523 


38.5174 


11.9112 


0.774374 



Table 3: Coefficients of the fall of FF for LIGO Livingston detector under the transforma- 
tion represented by equations (fT4l) and (fT31) for different 8 and f a . 
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Figure 4: Variation of FF with (f) T . 




5 10 15 20 



0.8 



0.6 



0.2 



a = a = n fw 


\\\ 

- w \ \ X 


e = e T = i° 

= T = 5° 

9 = S T = 10° 


■ \\ \ . 


= T = 15° 


\ 0-^0° 
\ T = 160° 


- \) \ x 

». 1 * 


\ f D = 1 Hz 


: \ \ \ 









5 10 15 20 



Figure 6: Fall of FF with T obs for different 6. 
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4 Conclusion 



In view of blind all sky search for CGW, we have studied the matching of signals from 
different source locations assuming the noise to be stationary and Gaussian. For fixed 
f , we observed that the matching of signals from locations in the sky that differ in their 
colatitude by n radians is independent of T obs , and 0. However, it falls with f Q . But in 
longitude the matching of signals falls with T obs , 9, (j), f Q . We believe that the matching of 
signals will increase for the real data. This is due to the resolution provided by the Fast 
Fourier transform (FFT). However, it may not be relevant to account this symmetries for 
the search of CGW from the output of LIGO-I, GEO600 and TAMA detectors. 

This analysis will be more relevant if one performs hierarchical search (Mohanty & 
Dhurandhar 1996; Mohanty 1998). This search is basically a two step search, in first step 
the detection threshold is kept low and in the second step a higher threshold is used. The 
higher threshold is used for those templates which exceed the first step threshold. The 
study of the matching of signals in different source locations will be also relevant for Laser 
interferometer space antenna (LISA) (Hough J., 1995). The work has been initiated and 
may be useful for the data analysis of CGW. 
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